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Abstract. 

Galactic short period double white dwarfs (DWD) are guaranteed gravitational wave (GW) 
sources for the next generation of space-based interferometers sensitive to low-frequency GWs 
(10 -4 - 1Hz). Here we investigate the possibility of constraining the white dwarf (WD) proper- 
ties through measurements of apsidal precession in eccentric binaries. We analyze the general 
relativistic (GR), tidal, and rotational contributions to apsidal precession by using detailed He 
WD models. We find that apsidal precession can lead to a detectable shift in the emitted GW 
signal, the effect being stronger (weaker) for binaries hosting hot (cool) WDs. We find that in 
hot (cool) DWDs tides dominate the precession at orbital frequencies above ~0.01 mHz (~1 
mHz). Analyzing the apsidal precession of these sources only accounting for GR would poten- 
tially lead to an extreme overestimate of the component masses. Finally, we derive a relation 
that ties the radius and apsidal precession constant of cool WD components to their masses, 
therefore allowing tides to be used as an additional mass measurement tool. 
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1. Introduction 

The Laser Interferometer Space Antenna (LISA) is expected to provide the largest ob- 
servational sample of short period DWDs (jNelemans et al. 20011 2001b, 2004. ILiu et al. 20101 
Rui ter et al. 2010)) . Recent theoretical calculations by Willems et al. (2007) and Thomp- 
son (2010) predicted a population of eccentric DWDs. Willems et al. (2007) studied eccen- 
tric binaries formed through dynamical interactions in globular clusters, while Thompson 
(2010) studied eccentric DWDs as products of hierarchical triple systems, where the ter- 
tiary induces Kozai oscillations in the inner binary (see also IGould 20111 for a discussion 
of the implications of such a population). Our focus here lies on these eccentric bina- 
ries, where measurements of apsidal precession could be used to constrain some of the 
WD properties. For this purpose, we use detailed He WD models with masses between 
0.1-0.3 Mq ( |Deloye et al. 2007] ) to investigate apsidal precession in eccentric binaries. In 
what follows, we consider a binary hosting two stars with masses Mi*t, and radii i?i,2, 
uniformly rotating with angular velocities fii^- We take the axes of rotation to be or- 
thogonal to the orbital plane. Let 7 be the argument of the periastron, P and a the 
orbital period and semi-major axis, respectively, and e the orbital eccentricity. We con- 
sider orbital and rotational periods long compared to the free harmonic periods of the 
component stars ( |Cowling 19 38, Sterne 1939, Sme yers &: Willems 2001[ ), so that we can 
neglect resonances between the tidal forcing angular frequencies and the eigenfrequencies 
of WD free oscillation modes, that is, we work in the static-tide limit. 
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2. Apsidal Precession and Driving Mechanisms 

Apsidal precession occurs when the binary's gravitational potential is perturbed from 
a pure Keplerian potential. This perturbation can be due, for example, to tides, rotation, 
and GR. The contribution to the apsidal precession due to the quadrupole tides raised 
in star i by the companion is (|Sterne 1939|) : 

307T (Rj\ 5 M 3 -, l + je 2 + |e\ 
W = -jT {-) — h ( i_ e2) 5 k ^ 

where ki is the so-called apsidal precession constant for star i, and it measures the star's 
central concentration. We refer to Valsecchi et al. (2011) for details on how to calculate 
ki from detailed WD models. 

The apsidal precession due to the rotational distortion of the star caused by centrifugal 
forces takes the form (e.g. Sterne 1939): 

2ir fRA 5 Mi + M 2 (Qi/O) 2 



A r**,i = -p{-) Mi (1 _ e 2 )2 **. (2-2) 

where Vt = 2tt/P is the mean motion. Finally, the GR contribution to apsidal precession 
at the leading quadrupole order is given by (jLevi-Civita 1937|) : 

/2tA 5/3 3G 2 / 3 (M 1+ M 2 ) 2 / 3 
laR=[ T ) — (1 _ e2) , (2-3) 

where G and c are the gravitational constant and the speed of light, respectively. The 
total apsidal precession rate (7) is the sum of the three contributions above considering 
both components. We note that all contributions to 7 depend on P, e and Mj (with jgr 
being dependent upon the total system mass), while only 7R t,j and ^Tid,i depend on Ri 
and ki. 
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Figure 1. Apsidal precession rates for cool (left) and hot (right) DWDs with masses 
Mi = M2 = 0.3 Mq, and orbital eccentricity e =0.01. The tidal and rotational rates include the 
contribution of both stars. The rotational contribution is calculated assuming synchronization 
at periastron. For cool WDs jridj and jR t,% decrease because of their dependence on iif . 

In Fig. [1] we show an example of tidal, rotational, and GR apsidal precession rates 
calculated with our models, as a function of the orbital frequency v = l/P. In gen- 
eral, GR dominates 7 at small v, while tides and rotation become more significant as v 
increases. The currently observed sample of rotation rates in DWDs suggests that ap- 
sidal precession is tidally-induced at frequencies where GR is no longer important (see 
Valsecc hi et al. 20111 for details). Apsidal precession can be significant for both hot and 
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cool DWDs. For cool components at v ~ 1 mHz, and e = 0.01, apsidal precession is ex- 
pected to shift the GW signal by 'Sn-An over 1 yr mission, which is detectable in proposed 
space based interferometers fsee I Willems et al. 20081 for a detectability study). The shift 
increases by more than a factor of 10 for hot WDs, and at v ~ 1 mHz tides dominate. 
Note that, due to the dependency of 7rid,i on Pf , binaries hosting hot (and therefore 
bigger) WDs have 7™ systematically higher than cool (and contracted) WDs by a factor 
of - 20. 

With lower component masses, "/gr,i decreases, while jridj and 7Rot,i increase. The 
larger radii of the lower mass models act to increase the tidal and rotational contributions 
relative to GR, therefore making detection even more favorable. Note that theoretical cal- 
culations from Willems et al. (2007) predicted orbital eccentricity distributions (affected 
by post-formation GW emission) that typically peak at e ~0.1, with a non- negligible 
fraction of systems having e up to ~0.7. Increasing e would lead to much higher rates of 
apsidal precession (see IValsecchi et al. 2011|) . 

4. Tidal Contributions to Apsidal Precession 

Given the dependence of "fen only on P, e, and Mi, if apsidal precession is detected in 
binaries where GR is the leading mechanism, it can be used to place constraints on the 
total system mass. Here we show that jTid,i can be used in a similar fashion considering 
the evolution of the term fc^P^ entering Eq. (|2.1[) as the WD cools. As Figure [5] shows, 
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Figure 2. Evolution of fci-Rf as a function of the WD's age. Filled squares denote plateau 
values. The maximum age considered for each model is 10 Gyr. 

during most of the WD's lifetime A^Pf assumes a plateau value which is unique for each 
mass. Consequently Eq. (I2.1[) effectively depends only on P, e, and Mi for much of the 
WD lifetime. A fit of the plateau values of fc^Pf as a function of the WD mass yields 
kiRf = -0.632 + 0.370 • M^ 1 '™, where M l is the mass in solar units, and fcjPf is in 
10~ 10 Pq . Therefore, for any DWD for which apsidal precession is detectable, 7 yields 
constraints on a particular combination of the masses in the system, which we call the 
" apsidal mass function" . 

5. Misclassification of the Sources 

Analyzing an apsidal precession detection only accounting for the GR contribution, 
could potentially lead to a dramatic overestimate of the components masses. As shown 
in Figure [21 an eccentric DWD could be misclassified as a binary black hole or neutron 
star if the total apsidal precession is mis-attributed to GR alone. 
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Figure 3. Lines of constant apsidal precession rate calculated for cool DWDs with v — 1 mHz, 
and e = 0.01. The left panel depicts lines of constant total j calculated for different combinations 
of components' masses accounting for GR, rotation and tides. The right plot shows, for these 
same 7, the masses that would be inferred solely accounting for GR. The various lines correspond 
to 7 in deg/yr of ~ 22570 (black-solid), ~ 5700 (dotted), ~ 2200 (dashed), ~ 1100 (dot-dashed), 
~ 680 (grey-solid). 

6. Summary and Conclusions 

We use detailed He WD models to investigate apsidal precession due to tides, rotation 
and GR in eccentric DWDs. We find that apsidal motion can lead to a significant shift in 
the emitted GW signal, the effect being stronger for binaries with hot WDs. In general, 
GR dominates the precession at small orbital frequencies v, while 7 is driven by tides 
at larger frequencies, where GR is no longer important. We investigate the astrophysical 
information that can be extracted from a measured 7. We show that, as the WD ages, 
fci-Rf reaches a plateau value that depends only on the WD mass. This behavior allows us 
to simplify the form of Eq. (|2.ip . and to use "jTidj to place constraints on the apsidal mass 
function also for short-period binaries, extending the possibility of mass measurement to 
the entire LISA frequency domain. We also find that the component masses that would be 
inferred from the GR contribution alone could be overestimated by orders of magnitude. 
We conclude that the inclusion of tides is necessary for the proper analysis of detected 
apsidal motion rates of DWDs in the LISA band. 
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